NASA CONTRACTOR 
REPORT 


TT 


CN 

CNI 


I 



<c 


<c 



AN ANALYTICAL STUDY OF THE RESPONSE 
OF A CONSTANT-ATTITUDE AIRCRAFT 
TO ATMOSPHERIC TURBULENCE 


by Frederick 0. Smetana and Ronald K. Carden 


Prepared by 

NORTH CAROLINA STATE UNIVERSITY 

Raleigh, N.C. 27607 

for Langley Research Center 


CR-2204 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • MARCH 1973 



1. Report No. 2. Government Accession No. 

NASA CR-2204 


4. Title and Subtitle AN ANALYTICAL STUDY OF THE RESPONSE 
OF A CONSTANT-ATTITUDE AIRCRAFT TO ATMOSPHERIC 
TURBULENCE 


7. Author(s) 

Frederick 0. Smetana and Ronald K. Carden 


9. Performing Organization Name and Address 


North Carolina State University 
Raleigh, N.C. 27607 


-■} j yUSTKi 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. Supplementary Notes 


3. Recipient's Catalog No. 


5. Report Date 

March 1973 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 

736-05-10-02 


11. Contract or Grant No. 

NAS I -9603 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 


16. Abstract 

A light airplane equipped with an automatic control system which drives large wing 
flaps and the stabilator so as to produce a constant pitch attitude in all flight modes is 
analyzed for its response to a specific gust. The aircraft is also equipped with a bank-angle 
steering, zero sideslip automatic control system which is studied for its effectiveness in 
suppressing a specific lateral gust. The gusts are assumed to be comprised of 200 lateral 
and 400 vertical sinusoids. Each is used to excite the controlled aircraft and the time 
response to the sum of all sinusoids is plotted. The assumption is that the gust may be 
treated as stationary in space but variable in time rather than the reverse as is usually the 
case. Results indicate that such a control system can suppress vertical gusts up to the limit 
of control authority. Either the lateral accelerations or the yawing velocity response to 
lateral gusts can be suppressed with this system but not both simultaneously. 


17. Key Words (Suggested by Author(s) ) 

Gust response 
Automatic control system 
Light aircraft 


19. Security Classif. (of this report) 

Unclassified 



20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

50 


22. Price* 

$3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 

























Page Intentionally Left Blank 



TABLE OF CONTENTS 


LIST OF TABLES 

LIST OF FIGURES 

NOMENCLATURE 

I NTRODUCT ION 

METHOD 

AIRCRAFT TRANSFER FUNCTIONS 
Latera l-D i rectiona I 

Longitudinal 

DISCUSSION 

CONCLUSIONS 

LIST OF REFERENCES . . . . 



LIST OF TABLES 


Page 


Table I. Fourier coefficients for lateral gust 20 

Table II. Fourier coefficients for longitudinal gust .... 22 

Table III. Open loop lateral transfer functions 26 

Table IV. Lateral equations of motion parameters 27 

Table V. Longitudinal equations of motion parameters .... 28 

Table VI. Open loop longitudinal transfer functions 29 


i v 



LIST OF FIGURES 


Page 


Figure I. Time history of vertical gust used to excite 

aircraft 30 

Figure 2. Time history of lateral gust used to excite 

aircraft . . 31 

Figure 3. Approximate one-dimensional power spectral 

density of vertical gust 32 

Figure 4. Approximate one-dimensional power spectral 

density of lateral gust 33 

Figure 5A. Open loop angle of attack response to vertical 

gust input 34 

Figure 5B. Closed loop angle of attack response to 

vertical gust input 35 

Figure 5C. Open loop pitch angle response to vertical 

gust input 36 

Figure 5D. Closed loop pitch angle response to vertical 

gust input 37 

Figure 6A. Open loop aerodynamic sideslip angle response 

to lateral gust input . 38 

Figure 6.B. Closed loop aerodynamic sideslip angle response 

to lateral gust input . . 39 

Figure 6C. Open loop roll angle response to lateral gust 

i nput . . 40 

Figure 6D. Closed loop roll angle response to lateral 

gust i nput 41 

Figure 7. Unity feedback representat ion of lateral 

control system 42 

Figure 8. Modified aircraft ... 43 

Figure 9. Unity feedback representation of longitudinal 

control system 44 


v 



NOMENCLATURE 


Fourier Cosine Coefficients 
Fourier Sine Coefficients 
variable in Laplace domain 
period of fundamental frequency 
aircraft initial forward speed 
forward speed 

lateral gust velocity input 
sides! ip velocity 
longitudinal gust velocity input 
a loop feedback function 
3 loop feedback element 
<j) loop feedback element 
0 loop feedback function 
angle of attack 

vertical gust angle of attack input 

total angle of attack 

sides I ip velocity 

s i degust ang le i nput 

tota I s i des I i p ang I e 

general input 

aileron input 

flap input 

rudder input 

stabilator input 



<p 

ro 1 1 ang 1 e 

b 

to+a 1 roll ang 1 e 

e 

pitch angle 

0 T 

total pitch angle 


yaw angle 

to 

frequency 


compensation, an artificial modification of a control surface deflection 
command such that the apparent vehicle dynamics are altered 
a device in a control system by which one can remove — or 
compensate for — undesired dynamic character i st ics in the 
vehicle. The command may be modified by filtering and/or 
by adding signal components proportional to vehicle 
position, velocity, and acceleration to the command. 

Simple feedback of position is not considered compensation. 


shaping network, a device for automatically altering a command signal 

such that some aspects of a vehicle’s dynamic character- 
istics will be more heavily excited than others. 
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INTRODUCTION 


A recent design study (Ref. I) took as its problem the development 
of a simple-to-f ly, constant-att i tude light aircraft. The concept which 
evolved employs fast-acting, full-span, Fowler-type wing flaps — operated 
together for lift modulation at all flight conditions and differentially 
for roll control — a completely separate fuselage leveler driving the 
elevator, and a turn coordinator controlling the rudder and differential 
flap motions. A flap-power interconnect is also provided to permit ascent 
at constant forward speed and airspeed changes at constant rate-of-c I imb. 
The pilot has only two controls: a wheel and a throttle. Changes in a 

throttle position change speed only. Sustained rotation of the wheel 
produces a coordinated turn. Forward motion of the wheel causes descent 
at constant forward speed and rearward positioning of the wheel results 
in a rate-of-c I imb at constant forward speed. 

The effort to provide the aircraft with satisfactory handling 
characterist ics led to the specification of powerful, fast-acting servo- 
controls. Such controls, however, can also form the basis of a gust 
alleviation system. With safety of flight in turbulent weather an area 
of growing national concern it was felt that an investigation of the 
applicability of the constant attitude concept to gust alleviation ser- 
vice was indicated. 

The availability of a digital computer program (Ref. 2) to calculate 
and plot aircraft time histories given involved transfer functions made 
possible its frequent use in visualizing the response of an aircraft to a 
variety of control imputs. The ease with which such results are inter- 
preted when contrasted to the transfer function representation led to a 
search for ways in which analytical gust responses could be presented in 
the time domain. (It is, of course, conventional practice to represent 
gust response as the product of the appropriate gust transfer function 
amplitude ratio squared and the power spectral density of the gust.) This 
approach would also facilitate direct comparison with flight records. 

It is generally assumed that the air motions associated with a gust 
are stationary in time but vary from point to point in space. Thus the 
aircraft will encounter angles of attack and sideslip due to gusts which 
vary across the dimensions of the aircraft and with changes in flight ve- 
locity. One common simplification to this picture is to assume that a 
vertical gust, for example, does not vary in the span-wise direction. It 
is therefore a one-dimensional gust. This gust is the one most frequently 
used in analysis because of its simplicity. It has been stated (Ref. 3) 
that the effects of encountering different portions of the gust wave with 
different portions of the aircraft — so called penetration effects — can 
usually be treated adequately by including the effect in the stability de- 
rivatives involving a or w. If this is the case, it would seem reasonable 
to conclude that one could obtain similar results by assuming the gust to 
be stationary in space and variable in time. One would merely need to 



include a time delay term — proportional to flight velocity — in these deriv- 
atives. The stationary-in-space, variable-in-time model could be used in 
any event for gust wave lengths much larger than the distance between the 
wing quarter-chord and the tail quarter-chord . For most light aircraft 
a gust frequency of less than 1.0 to 1.5 radians/sec will insure that this 
condition is satisfied. 

The time response of a linear, uncoupled aircraft to a sinusoidal 
disturbance of a control surface is usually fairly easy to determine. 

Since the system is linear, response to complex control surface inputs can 
be determined by superposition. Then if a gust can be considered to 
represent a disturbance similar to that induced by a control surface de- 
flection, the time response of the aircraft to a gust input can be com- 
puted in a similar fashion. This is the viewpoint adopted for the present 
work. 
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METHOD 


Flight records of aircraft response to gusts have the unfortunate 
aspect that the gust exciting that particular response is almost never a 
"standard" gust. The principal differences will be in mean square ampli- 
tude and in harmonic content at various frequencies. However, unless the 
aircraft has highly resonant stability modes, the departure in harmonic 
content with frequency from that of theoretical models or standards is a 
second order effect compared with the variation in amplitude so far as 
aircraft response is concerned. This will be shown aposteriori. Ampli- 
tude variations of course are easy to standardize so long as the system 
is linear. Thus it would seem that any we I i -documented sample of atmos- 
pheric turbulence would provide a suitable forcing function. 


A recent report (Ref. 4) contained measurements of vertical and 
lateral gust velocities taken at relatively high altitudes during a storm 
in the Rocky Mountain area. Since the gust velocities were rather severe 
the records seemed to be well-suited for determining the harmon i c. content 
of the gusts. Accordingly, the period of time between about 15 seconds 
and 90 seconds (see figures I and 2) was selected as the fundamental 
period. The curves were photographically enlarged and the gust amplitude 
read at 638 evenly-spaced time increments. The lateral gust amplitude was 
read at 572 evenly-spaced time increments. Each of these data sets was 
then expanded by a factor of 4 by interpolating linearly between suc- 
cessive data points. The 2549 data points comprising the vertical gust 
and the 2285 data points comprising the lateral gust were processed by a 
Fourier analysis routine to obtain cosine amplitude, A w , and sine ampli- 
tude, B , at each frequency. Approximately MOO harmonics of 2tt/75 were 
computed. f(t) was computed from the series and compared with the orig- 
inal photograph. It was found that 200 harmonics were sufficient to re- 
produce all except the very sharpest peaks of V* while 400 harmonics were 
necessary to obtain a good fit to W*. Table I presents A^, and for 
the 200 frequencies needed to form Vg. Table II presents the same re- 
sults for W* These tables give the amplitudes and frequencies of the 
sine and cosine forcing functions used as the imputs to the aircraft trans- 
fer functions. The quantity 


h (A o ) 2 + 


can be shown to approximate closely the one-dimensional power spectral 
density of the gust. Figures 3 and 4 show 



+ B (j 2 


) 


plotted against log frequency. Except for the scatter the plots are some- 
what similar to those for the Dryden turbulence model. 
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This model predicts that the power spectral density will fall off as 
frequency to the -2 power above a certain frequency. The present data 
show a slope of about -3. Note also that in the present case the verti- 
cal and lateral gust amplitudes and harmonic content are somewhat dif- 
ferent, indicating that the turbulence was not isotropic. 

Measurements made during flights in clear air turbulence over the 
Eastern Seaboard (Ref. 8) show a slope of about -2 at the higher fre- 
quencies. Scatter in these data ranged from about ±15$ at high fre- 
quencies to about ±40$ at low frequencies. The data in Ref. 8 covered a 
frequency spectrum from about half the minimum value used here to about 10 
times the maximum value used here. The low frequency amplitude of the 
gust analyzed here was about an order of magnitude larger than that found 
in Ref. 8. At the middle frequencies, n, | Hz, the amplitudes were approx- 
imately the same while at the higher frequencies the amplitudes in the 
present case were less than those reported in Ref. 8. The fact that the 
gust time histories of Ref. 4 were inferred from aircraft responses and 
plotted in such a way as to obscure very high frequency information is 
probably responsible for the more rapid attenuation in gust amplitude 
found at frequencies above 0.13 Hz by the Fourier analysis. 

It may be noted that for the present calculation the frequencies of 
the aircraft stability modes (except for the rolling mode) are all less 
than I Hz. Thus, this more rapid attenuation in gust amplitude with fre- 
quency has no significant effect on the displacement of the airframe cen- 
ter of gravity but rather will be felt primari ly as a somewhat reduced 
"shaking” of the airframe. 

The time response of the system to the function 

where go = — ( I ) 

T 


was found by multiplying the aircraft transfer functions by the Laplace 
transform of equation (I) and then solving for the time response by the 
method of residues. The time responses of the system for each value of 
n, n - 0, I ,2, . . .200, were then summed to find the total system response 
to a lateral gust. A similar procedure was used for the vertical gust. 

The results are shown in figures 5 and 6 . It may be seen that within 
the assumptions used to develop the aircraft transfer functions (see next 
section for development) the vertical gust alleviation produced by the 
constant attitude control system is significant. The side gust allevia- 
tion appears to be outstanding. The system does in fact reduce the 3 and 
( p errors to very small values. But because the air mass is moving rela- 
tive to the aircraft, the gust alleviation capability is actually somewhat 
more limited than it appears as will be explained in a later section. 
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AIRCRAFT TRANSFER FUNCTIONS 


Latera I -Pi reef iona I 

The analysis of the response of an aircraft with rudder and aileron 
control loops to a lateral gust input has been treated in a very general 
fashion by Onstott and Salmon (Ref. 5). In terms of the present problem, 
the output variables chosen for minimization are <j> and 3- In their treat- 
ment both 3 and (f> must be measured in an inertial reference system. This 
character i zat i on of 3 is unfortunate because for simplicity and reliability 
one would like to use a vane or other aerodynamic device to sense ft. A 
vane, of course, measures the orientation of the aircraft body axis rela- 
tive to the instantaneous air mass velocity over the aircraft and thus 
does not give ft in an inertial reference system. This is easily seen from 
the following example. 

Consider an aircraft with stable I atera I -d i rect iona I dynamics to be 
flying along a straight line at constant speed in still air. Suddenly, a 
step gust approaches from the right. Before the aircraft can respond, 
it is subject to a positive ftg which will be sensed by a vane, ft as 
given by Onstott and Salmon is zero because the aircraft has not yet 
turned with respect to the original heading. After a while the aircraft 
will relieve the side force imposed on it by the gust in two ways: by 

accelerating its c.g. to the gust velocity and/or by rotating so as to 
align itself with the resultant air mass velocity. When all the system 
dynamics have decayed, the vane will indicate zero sideslip, but the sum 
of the ground track angle and the heading wi II now correspond to the vec- 
tor sum of the aircraft velocity and the gust velocity. In the inertial 
reference frame there will be a steady state final value of ft equal to 
-ftg. This is what the equations of Onstott and Salmon predict. Thus if 
one wishes to use a vane as the feedback control element, he must rewrite- 
the equations to account for its behavior. 

The need for this rewriting makes it attractive to use a somewhat 
simpler derivation than that presented by Onstott and Salmon to obtain the 
required result. Since it is assumed that the system is linear, a sum of 
elementary solutions is also a solution. It is a familiar school exercise 
to obtain 4>, ft, and ip in response to either a rudder or aileron deflection. 
(See for example Ref. 6 for elaboration.) The response of a system to both 
inputs simultaneously is simply a sum of the response of the system to the 
inputs separately: 

% - <^> S A + <fj> «R <2> 
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= ( f^ 6 A + ( f^ 6 R 


(3) 


The quantities in parentheses represent the indicated open loop 
transfer functions. (See standard texts or Ref. 6 for development.) The 
various open loop transfer functions appear in Table III with numerical 
values substituted. <J> T and ftp are the total response to combined rudder 
and aileron inputs. The derivation of the transfer functions in Ref. 6 
assumes the air mass to be stationary. Hence, <j>j and By in equations (2) 
and (3) are measured in an inertia I reference system, even though B is 
usually thought of as an aerodynamic quantity. The response to a gust 
input may be included in the total response simply by adding the term 

‘t;’ 


to the bank angle equation and the term 

Bn 


( r" ) + ' 

LV J 


to the sideslip equation-. Since <p is always measured in an inertial 
reference system the open loop expression developed by Onstott and Salmon 
for 

( 4-, 

V 


may be used unaltered. However, in order to describe the value of By 
which would be sensed by a vane mounted on the aircraft, it is necessary 
to add the sidewise air mass velocity to the sidewise aircraft velocity 
in inertial space as computed by Onstott and Salmon; hence Bg is added 
to the right hand side of the equation for Bt in addition to The term 


(!-> B 


g* 


For the purposes of this analysis it is assumed the pilot makes no 
inputs to the control system during a gust encounter. The latera l-di rec- 
tional control system used for the study is shown in figure 7. The basic 
aircraft is shown in three views in figure 8 and its stability deriva- 
tives are given in Table IV. Because no pilot inputs are assumed the 
feedback elements i n the aileron and rudder loops can be represented by 


* <0- 1 


17) (S + 3) 


fjLi-') 

\S + 50/ 


(4) 
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Y s = (3.2HS ♦ 3.4XS 4 3.4) (rTToo) (5> 

respectively. Since the signals processed through these elements are fed 
into the aircraft negatively, one can represent the aileron and rudder in- 
puts to the aircraft as 

= ~ Y <j). h (6) 

6 R = _y 3 • (7) 

With these substitutions, the equations for <j>y and 3j may be solved simul- 
taneous I y to yield 



It can be shown that the denominators of these equations are identical to 
those one would obtain using the procedure of Onstott and Salmon. The 3j 
numerator, however, contains two additional terms, 

' i+ ^) ; 

and the <f>-|- numerator has the additional term 

‘ Y « (fe) 
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to account for the fact that an aerodynamic 3 is sensed rather than an 
i nert i a I . 

The responses 3y and <j>y are computed by evaluating 



by the method of residues, oj = and U Q is aircraft flight velocity. 

Responses for both open and closed loop cases are shown in figure 5. 

To accomplish the gust response attenuation noted in figure 5, the 
aircraft is yawed sufficiently to develop a side force equal and opposite 
to that imposed on the aircraft by the gust. In turbulent air where the 
gust approaches the aircraft first from one side and then the other in 
fairly rapid succession, fhe oscillations of the aircraft will be much 
larger than if the control system were not operational. However, the 
oscillations are phased so that the lateral acceleration of the aircraft 
c.g. is almost zero and ground track remains straight. 

If one wished instead to exchange oscillations about a c.g. moving 
along a straight line for lateral accelerations of the whole airplane and 
an altered ground track this is easily done by measuring inertial 3 in- 
stead of aerodynamic 3. The control system will then deflect the rudder 
only enough to develop sufficient yawing moment to keep the aircraft from 
rotating (but not sufficient to rotate the fuselage so that it develops a 
balancing side force.) While the heading remains essentially constant, 
the fuselage is accelerated first to one side and then to the other. Be- 
cause of the aircraft inertia, however, the fuselage will not move as 
rapidly as the impinging gust. 

The construction of an inertial 3 sensor is somewhat more complex 
than is that of a simple aerodynamic vane. From the equations of motion 
one sees that 
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( 12 ) 


I 

/ 


3 = 


_ o 


(a v - g<J>)dt 


U, 


Hence one requires a lateral accelerometer, a measurement of airspeed, 
and an indication of bank angle. The bank angle signal multiplied by a 
constant is subtracted from the accelerometer signal in a summing am- 
plifier. The resulting signal is then integrated and fed to one terminal 
of a divider. The other terminal is fed with a signal proportional to 
airspeed. The output of the divider is then the inertial 3 signal. 

There is no clear-cut human factors evidence to support a choice of 
which aircraft motion is inherently preferable in response to a gust. 
Probably the magnitudes of the two motions will be the determining fac- 
tor. If the sidewise displacement of the c.g. is fairly small it would 
seem to be preferable to the aggravated rotation necessary to neutralize 
the aerodynamic side force. During approach, however, one would probably 
prefer a constant ground track. 

It may be noted that the airframe dynamics are the same with either 
aerodynamic or inertial 3 input. The same compensation circuit may also 
be used. (Compensation is defined on page vili.) 

There is the problem, however, that the indication of an inertial 3 
sensor constructed in the manner described above will drift over long 
periods of time. To use this sensor successfully, it will probably be 
necessary to pass the output signal through a high pass filter and sum 
the result with the output of a low pass filter fed by the position in- 
dication of an aerodynamic vane. The filter cross-over frequency should 
be well below the natural frequency of any airframe oscillatory mode and 
in fact as low as the drift level permits. Rigorously, the use of a vane 
to provide D.C. error correction is possible of course only under the 
assumption that the gust has no components with frequencies below the 
cross-over frequency. Otherwise, the vane will force the aircraft to 
track these low-frequency components. 

There is also the possibility of sensing a directly rather than having 
to calculate inertial 3* The compensation required with acceleration in- 
put has not been investigated, however. 


Long i tud i na 1 

The longitudinal control system as designed modulates thrust so as to 
keep forward speed constant during ascent and descent unless throttle posi- 
tion is changed. For the purpose of this analysis, however, the thrust 
will be assumed to be constant. The flap motion is then controlled 
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entirely by the a sensor and stabilator motion by the 0 sensor. This 
simplified longitudinal control system is depicted in figure 9. In a 
manner analogous to the analysis of lateral motions 



The notation used in these equations follows that of references 2 and 6. 
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( 19 ) 


(S - X u )u - (SX^ + X w )a/U D - (SX q - g)0 = Xg^g + Xg^Sp + X a ^a g 

-(Z u )u + [SCI - Z-) - Z w ]a/U Q - S(U Q + Z q )6 = Zg 6 S + Zg 6p + Z a a g (20) 

o r 9 

-(M u )u - (SM^ + M w )a/U 0 + (S 2 - M q S)0 = Mg 65 + Mg 6 p + M a a g (21) 

S F g 


From these equations two general transfer functions may be derived: 


AaS 3 + BaS 2 + CaS + Da 

AgS 2 + BgS + C 0 
0 Dj 

where A a = Z 5 /U 0 

Ba = £-Z< 5 (M q + X u ) + MgU 0 ] /Uq 
C a = [x u (Z 6 M q - M 6 U 0 )] /U 0 
D a = [9 (z 6 M u - m 6 z u>] /Uq 

Aq = (ZgM^ + Mg) 

B 0 = Z 6 <M W - M^X U ) - Mg (X u + Z w ) 

c 6 = z 6 (M u X w - MwV + M 6 (x u z w - x w z u) 

D| = AS 1 * + BS 3 + CS 2 + DS + E 


( 22 ) 


( 23 ) 



with 


A = I 

B = - <M q + X u + Z w + U 0 M*) 

C = M q Z w - U 0 M W + X u (M q + Z w + U Q M*) - X W Z U 
D = - X u (M q Z w - U q M w ) - M u U q X w + M q X w Z u + g(M-Z u + M u > 
E = g(M w Z u - M U Z W ) 


The six individual transfer functions needed in the analysis are 
obtained from these two general functions as follows: 

(I) For- (~~) substitute X^, , etc., for Xg, Zg, etc., in the 

numerator of (~) 


(2) For (“r— ) substitute Xg , Z^ , etc., for Xg, Zg, etc 


S U S 


(3) For (2L.) substitute X^ , Z^ , etc., for X*, Z*, etc 

ut ^-1 CX o o 


g g 


A similar procedure is employed to find the 0 transfer functions. 

The values of the longitudinal stability derivatives, X^, Zg^_, 


etc., are given in table V. Table VI presents numerical values of the 

(Y 

transfer functions (-r— ) etc., as a ratio of polynomials. 

r 

The reader unfamiliar with this notation, the correspondance between 
the dimensional stability derivatives given here and the more widely used 
non-dimensional stability derivatives, and the equations of motion in the 
frequency domain rather than in the time domain is directed to Reference 6 
for a very lucid, comprehensive explanation. Those finding Reference 6 
difficult to locate may find NASA Contractor Report 1975 more accessible. 
This work treats many of the same questions in a similar fashion. 
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DISCUSSION 


It has been pointed out in the literature (Reference 3) that ef- 
fective gust suppression in the longitudinal mode requires control of both 
plunging and pitching motions. In the control system discussed herein 
pitching motion is restrained by stabilator operation and plunging motion 
by modulation of flap position. Unfortunately, in the lateral-di rectional 
mode there is no simple means to modulate the side force on the airframe 
without simultaneously yawing the fuselage. One can therefore suppress 
either the sidewise velocity or the yawing motion, but not both. Ex- 
perience with vertical gust alleviators has shown that suppression of only 
one of the two motions involved is not very satisfactory. It would appear, 
therefore, that completely satisfactory lateral gust suppression using 
this concept will require a substantial increase in vertical tail and 
rudder area, a decrease in effective dihedral (so as to develop a yawing 
moment due to rolling counter to the yawing moment developed by the rudder 
as it produces a side force to oppose the gust), and the use of 3 in- 
ertial, or perhaps ay, as the control system input. As now constituted, 
that is, with a vane 7 input the lateral- control system provides improved 
capability to fly a given ground track in exchange for somewhat greater 
osci I latory response to a side gust. 

Another possible means of lateral gust suppression which has not been 
investigated in detail is the use of combined rolling and pitching motions 
to develop the necessary side force. There are some indications that air- 
crew members might find such motions preferable to either yawing motions 
or lateral accelerations. A control system to provide the necessary mo- 
tions would have to couple the I atera I -d i rectiona I system with the forward 
speed-rate of climb system and would be useful only during gusts. For 
this reason its cost is likely to be fairly high and therefore its employ- 
ment on general aviation craft is less likely in the near future. 

A third possibility is the use of movable vertical surfaces mounted 
on the wings to modulate the side force on the aircraft. Two surfaces, 
each with a 2 r chord and a 3 ’ height, will provide the necessary side 
force generation capability for a typical 4-place general aviation craft. 
The surfaces would be actuated by the control system upon receipt of an 
(ay - g<f>) indication. The vertical tail would then be responsible only 
for control of yawing motions. An aircraft employing side force surfaces 
is discussed in reference 7. 

The suppression of vertical gust responses as shown in figure 5 is, 
in some respects, very satisfactory and in others less than one would wish. 
Note that the pitching response has been completely eliminated while the 
vertical motion has been reduced about 50$. This can be explained by the 
fact that the pitch loop (fuselage leveler) has a relatively high gain and 
the stabilator a very high aerodynamic effectiveness. The plunge channel, 
however, has the aerodynamic ability to compensate for changes in angle of 
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attack of only about ± 10°. Now, when a vertical gust of 30 ft/sec en- 
counters an airplane moving at 180 ft/sec it is equivalent, so far as the 
airplane is concerned, to a 9.4 degree change in angle of attack. If the 
vertical gust velocity is greater or the aircraft speed slower, the air- 
craft has insufficient flap effectiveness to alter the wing lift enough to 
compensate for the change in lift due to change in angle of attack. The 
aircraft will then seek to plunge (up or down) to relieve the excess lift 
force . 

In the example used, the vertical gust velocities reached 120 ft/sec 
which at an aircraft speed of 200 ft/sec, corresponds to a change in angle 
of attack of 31°. As the aircraft accelerates in the direction of the 
gust velocity, the effective forcing function is reduced so that after a 
while the vertical velocity of the aircraft equals that of the gust. For 
this very severe example it is seen that the aircraft's inertia tends to 
limit the vertical velocity in response to a gust while the control system 
will negate about 10° of the induced angle. The remaining vertical ve- 
locity divided by the forward flight speed is the response quantity plotted 
against time (expressed as degrees) in figure 3. Because of the need to 
scale gains for commands up to 10° and because of the presence of a well 
damped high frequency oscillation at the gain level chosen, it is not 
practical to choose a higher gain for the control system. Thus, the level 
of gust attenuation shown represents about the practical limit for this 
system. 

The dominant time constant in the control system is about 0.6 secs. 
Since the gust components begin to diminish in strength as their periods 
fall below I second, it is evident that most of the gust energy can be 
absorbed by this system if the gust amplitude is not too large. While a 
somewhat shorter time constant (say 0.1 secs) would be desirable, there 
is not a simple way this can be achieved because of the location of the 
open loop zeros in the a/<5p transfer function, the system gain limits, and 
the presence of other poles and zeros along the negative real axis to the 
left of - 5.0. Further study of this problem would be desirable. 

A word about the possible mechanization of the compensation circuits 
is probably in order since those familiar with control systems will prob- 
ably recognize some unconventional practices here. The longitudinal con- 
trol system involves two compensators: 


0 : 


(S + 7) (S + 8) 



(S + 3) (S + 4) 
S 
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The second of these presents few problems if an acce I erometer is 
used to measure w. It is seen that in addition to w itself, a first and 
second integral are required. Since angular acce I erometers are avail- 
able as well as rate gyros and position gyros, the function in the 0 loop 
(S + 7 ) ( S + 8) is easily mechanized. (S + 5.8)/(S + 0.58) can be built 
entirely with passive elements. 

The lateral control system also contains two compensators: 

<f>: (S + 3) 

3: (S + 3.4) (S + 3.4) 


the first of these requires only a rate gyro and a position gyro. No 
instrument presently exists which can measure v. Thus it will be 
necessary to develop this signal electrically either by differentiating an 
acce I erometer output or by multiplying the signal (S + 3.4) by itself. 

It is felt that recent advances in solid state device technology will per- 
mit this to be done successfully. 

Finally, it may be mentioned that because the pitch attitude with 
the control system operating is always near zero, the vertical acceleration, 
a guantity of interest when evaluating vehicle riding qualities, is just 

UoOt. Examination of Figure 5 shows that while the control system does not 
reduce the peak accelerations a great deal it is effective in suppressing 
the low frequency components of the vertical acceleration. As a result, 
the acceleration exceeds a given value for a much shorter period of time 
with the control system on than with it off. 
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CONCLUSIONS 


1. The control system for the s i mp I e-to-f I y , constant-attitude air- 
craft can provide significant alleviation of both vertical and lateral 
gusts, provided aircraft velocity is sensed inertial ly. Such a control 
system seems capable of very beneficial application to low wing loading 
STOL machines for gust alleviation. 

2. The vertical gust alleviation is limited by the aerodynamic 
capability of the flaps to modulate aircraft C[_ and by the difficulty one 
has in compensating electrically for inherently slow airframe dynamics. 
Because of a very aerodynamica I ly effective stabilator and the use of 
high gain in the feedback and servo systems, pitching due to gusts is 
very heavily suppressed. 

3. Because of the inability of the aircraft and its control system 
to generate significant sideforce without yawing, one cannot suppress 
both sidewise motion and yawing simultaneously. Either can be suppressed 
separately, however, the choice being determined by which of the motions 
is less objectionable. 

4. The suppression of linear motions requires inertial sensing of 
aircraft velocity. These signals then surplant those of aerodynamic 
sensors in the feedback loops. 

5. The representat ion of a gust by a specific series of sine and 
cosine disturbances in angle of attack (or angle of sideslip) is an 
effective means of developing a suitable forcing function with which to 
excite the aircraft transfer function and thereby determine the time 
response of the aircraft to a specific gust. The method is easily modi- 
fied to treat (a) lags in the build up of aerodynamic forces and (b) 
"penetration" effects. By accepting additional complexity the method 
can be extended to permit spatial variations in gust amplitude. 
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TABLE I, Fourier Coefficients for Lateral Gust 
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TABLE II, Fourier Coefficients for Longitudinal Gust 
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13.7421 

0.0407 

- 0. 1395 

9.5856 

0.0626 

- 0.0466 

13.8270 

- 0.0918 

0. 1905 

9.6704 

- 0.3849 

- 0.0926 

13.9118 

0.0118 

- 0. 1789 

9.7552 

0.3031 

- 0.1759 

13.9966 

0.0813 

- 0.2117 

9.8400 

- 0.4082 

- 0.0478 

14.0814 

0.0550 

0.0432 

9.9249 

0.2622 

0. 1276 

14.1663 

- 0.0017 

- 0.0837 

10.0097 

- 0.1007 

- 0.1678 

14.2511 

0.2274 

0 . 1381 

10.0945 

- 0.2421 

- 0.2754 

14.3359 

- 0.0157 

- 0.0205 

10.1793 

0. 3957 

- 0.0321 

14.4207 

0.0218 

0.0738 

10.2642 

0. 1676 

0.0149 

14.5056 

0.0446 

- 0.0873 

10.3490 

- 0.2029 

- 0. 1435 

14.5904 

- 0.0609 

- 0. 0606 

10.4338 

0. 3094 

- 0.0576 

14.6752 

0.1477 

- 0.0132 

10. 5187 

0. 1199 

- 0. 1139 

14.7601 

- 0.0275 

0.0235 

10.6035 

0.0637 

0. 3280 

14.8449 

0. 1030 

- 0.0955 

10.6883 

- 0.3474 

0.0117 

14.9297 

0.0705 

0. 1042 

10.7731 

0. 1329 

- 0.3025 

15.0145 

- 0.0621 

- 0.0931 

10.8580 

0. 2299 

- 0.1611 

15.0994 

0.0846 

- 0.0399 

10.9428 

0.0631 

- 0.0489 

15. 1842 

0.0954 

- 0.0933 

11.0276 

0.2003 

- 0. 1269 

15.2690 

- 0.0309 

0.0503 

11.1125 

- 0 . 05 C 1 

- 0.4010 

15.3539 

0.1433 

- 0. 1737 

11.1973 

0.0042 

0.2810 

15.4387 

0.0878 

- 0.0046 

11.2821 

0. 1289 

- 0.1605 

15.5235 

0.1538 

0. 0040 

11.3669 

- 0.0272 

- 0.0294 

15.6083 

- 0.1142 

- 0.0001 

11.4518 

0.0503 

0.0013 

15.6932 

- 0.0959 

0.0745 

11.5366 

0.0485 

- 0 . 1015 

15.7780 

- 0.0028 

- 0.0789 

11.6214 

0. 1251 

- 0.0627 

15.8628 

- 0.0531 

- 0.0710 

11.7063 

0.0241 

0.1840 

15.9476 

- 0.0264 

0.0103 

11.7911 

- 0. 1481 

- 0.0284 

16.0325 

0. 1068 

- 0 . 1576 

11.8759 

0.0189 

- 0.0813 

16.1173 

0.0782 

0. 0673 

11.9607 

0. 1391 

- 0.2560 

16.2021 

- 0.0211 

- 0 . 1823 

12.0456 

- 0.0722 

0.3051 

16.28 70 

0.1134 

0. 1241 

12. 1304 

0.0254 

- 0.4272 

16.3718 

0.0156 

- 0.0558 

12.2152 

0. 1526 

0.0525 

16.4566 

0.2049 

- 0.0202 

12.3000 

0.2830 

0.0549 

16.5414 

- 0.1290 

- 0.0435 

12.3849 

0.1176 

0.0311 

16 . 6263 

0.1564 

- 0.0332 

12.4697 

- 0 . 1680 

0.0391 

16.7111 

- 0.0211 

0.0385 

12.5545 

- 0.0082 

- 0.0948 

16.7959 

0.0557 

0.0117 

12.6394 

0.0645 

- 0. 1480 

16.8808 

0.0086 

- 0. 0474 


TABLE II. (continued) 



(Jj 

A 

B 

0) 

A 



fa ) 

a) 


(u 

a* 

16.9656 

- 0.0331 

- 0.0804 

21.2070 

0. 1544 

- 0.0266 

17.0504 

- 0.0053 

3.0126 

21.2918 

- 0.0934 

0.0003 

17.1352 

- 0.0697 

- 0.0868 ; 

21.3766 

0.0003 

- 0.0595 

1 7.22 OL 

0.1032 

- 0. 1626 

21.4615 

0.0900 

- 0.0586 

17.3049 

0.0505 

0.0808 

21.5463 

0.0322 

0. 0269 

1 7.3897 

0. 1491 

- 0.0231 

21.6311 

0.0374 

0.0210 

17.4745 

0.0543 

- 0.0073 

21.7159 

0.0282 

- 0.0231 

17.5594 

0.0660 

- 0.0546 

21.8008 

0. 0054 

0.0581 

17.6442 

- 0.0756 

0.0051 

21.8356 

0.0963 

0.0999 

17.7290 

0. 1347 

0.0365 

21.9704 

- 0.0791 

- 0.0374 

17.8139 

- 0. 0602 

- 0.0069 

22.0551 

0.0266 

- 0.0651 

17.8987 

- 0. 0394 

- 0.0534 

22.1401 

0.0133 

( 7.0359 

17.9835 

0.0875 

- 0.0647 

22.2249 

- 0. 0387 

0.0732 

13.0683 

0.0169 

0.0792 

22.3097 

0.0795 

- 0.0881 

18. 1532 

- 0.0814 

0.0045 

22.3946 

- 0.0061 

- 0.0045 

18.2380 

0.0760 

- 0.0580 

22.4794 

0.0568 

0.0732 

18.3228 

0.0732 

0.0082 

22.5642 

0.0414 

0.0098 

16.4077 

0.0153 

- 0. 0184 

22.6490 

0.0054 

0.0005 

10.4925 

0. 1527 

0.0143 

22.7339 

0.0099 

- 0.0348 

18.5773 

- 0.0227 

- 0.0503 

22.8187 

- 0.0174 

0.0315 

18.6621 

0.0779 

0.0017 

22.9035 

0.0845 

- 0.0344 

18.7470 

0.0028 

0.0714 

22.9864 

0.0900 

- 0.0052 

16.8310 

0.0914 

0.0104 

23.0732 

- 0.0845 

- 0.0408 

18.9166 

- 0.0259 

- 0.0466 

23.1580 

0.0242 

- 0.0150 

19.0014 

- 0.0129 

- 0.0049 

23.2428 

0.0198 

- 0.0670 

19.0863 

- 0.0935 

- 0.1278 

23.3277 

0.0811 

- 0.0219 

19.1711 

0.1190 

- 0.0137 

23.4125 

- 0.0506 

- 0.0662 

19.2559 

- 0.0310 

- 0.0433 

23.4973 

0.0587 

- 0.0018 

19.3408 

0.0859 

0.0415 

23.5822 

0.0617 

- 0.0109 

19.4256 

- 0.0378 

- 0.0624 

23.6670 

0.0263 

- 0.0327 

19.5104 

0.1110 

0.0204 

23.7518 

- 0.0074 

0.0035 

19. 5952 

0.0161 

- 0.0825 

23.8366 

- 0.0002 

0.0048 

19.6801 

0.0696 

0.0080 

23.9215 

0.0648 

0.0945 

19. 7649 

- 0.0245 

0.0187 

24.0063 

- 0.0091 

- 0. 1039 

19. 8497 

0.0878 

0.0153 

24.0911 

0.0383 

- 0. 1079 

19.9346 

- 0.0223 

- 0. 0540 

24.1769 

0.0281 

0.0742 

20.0194 

- 0.0102 

- 0.0615 

24 .2608 

0.0578 

- 0.0140 

20.1042 

0.0263 

0.0576 

24.3456 

0.0423 

- 0.0254 

20.1890 

- 0.0915 

- 0.0712 

24.4304 

0.0253 

0.0312 

20.2739 

0.1294 

- 0.0518 

24.5153 

- 0.0154 

0.0165 

20.3587 

- 0.0234 

- 0.0308 

24.6001 

0.0947 

0.0294 

20.4435 

0.0447 

- 0.0384 

24.6849 

0.0472 

- 0.0513 

20.5284 

0.1147 

0.0453 

24.7697 

- 0.0017 

0.0834 

20.6132 

0.0913 

- 0.0304 

24.8546 

0.0149 

0.0729 

20.6980 

0.0143 

0.0161 

24 . 93^4 

0.0310 

- 0.0566 

20.7828 

- 0.0272 

0.0589 

25.0242 

o.oeo3 

C .0331 

20.867 7 

0.1353 

- 0.0813 

25.1091 

- 0.0449 

- 0.0374 

20.9525 

0.0824 

0. 0436 

25.1939 

- 0.0541 

0.0305 

21.0373 

- 0.0234 

- 0.0423 

25.2787 

0.0798 

- 0.0655 

21.1221 

- 0.0404 

0. 0045 

25.3635 

0.0938 

- 0.0253 


TABLE II. (continued) 
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(Jj 

A 

B 


<0 

U) 

25.4484 

- 0.0307 

- 0.0006 

25.5332 

0.0212 

0.0295 

25.6180 

0.0319 

- 0.0197 

25.7029 

0.0209 

0.0268 

25.7877 

- 0.0322 

- 0.0488 

25.8725 

- 0.0249 

0 . 1244 

25.9573 

0.0218 

0.0092 

26.0422 

- 0.0133 

- 0.0680 

26.1270 

0.0129 

- 0.0366 

26.2118 

0.0171 

0. 0788 

26.2966 

0.0160 

- 0.0027 

26.3815 

0. 0797 

- 0. 0740 

26.4663 

0.0241 

- 0.0332 

26.5511 

0.0887 

0.0708 

26.6360 

0.0548 

- 0.0070 

26.7208 

- 0.0095 

- 0.0343 

26.8056 

0.0291 

0.0281 

26 . 89 C 4 

0.0148 

0.0408 

26.9753 

0.0612 

0.0013 

27.0601 

0.0253 

- 0.0565 

27.1449 

- 0.0240 

- 0.0311 

27.2298 

0.0134 

- 0.0200 

27.3146 

0.0312 

0.0340 

27. 3994 

0.0464 

- 0.0429 

27.4842 

0.0156 

- 0.0791 

27.5691 

0.0454 

0.0445 

27.6539 

0.0776 

0.0122 

27.7387 

- 0.0260 

0.0011 

27.8235 

- 0.0449 

- 0.0340 

27.9084 

0.0346 

0.0481 

27.9932 

0.0178 

0.0510 

28.0780 

- 0.0174 

- 0.0285 

28. 1629 

- 0.0813 

- 0 . 039 / 

28.2477 

0.0392 

0.0614 

28.3325 

0.0149 

- 0.0294 

28.4173 

0.0406 

0.0013 

28.5022 

0.0022 

- 0.0156 

28.5870 

0.0081 

0.0259 

28.6718 

0.0828 

- 0. 0299 

28.7567 

0.0555 

0.0119 

28.8415 

- 0.0620 

0.0245 

28.9263 

0.0108 

0.0312 

29.0111 

0.0344 

- 0.0669 

29.0960 

0.0351 

0.0168 

29. 1808 

- 0.0382 

- 0.0459 

29.2656 

- 0.0184 

- 0.0237 

29.3504 

0.0715 

0.0141 

29.4353 

0.0279 

- 0.0447 

29.5201 

- 0. 0138 

0.0236 

29.6049 

0.0151 

- 0.0099 


TABLE II, (continued) 


OJ 

A 

B 


Xl ) 

t 0 

29.6898 

0.0460 

0.0106 

29.7746 

0.0330 

- 0.0046 

29.8594 

- 0.0109 

- 0.0156 

29.9442 

0.0762 

0.0621 

30.0291 

0.0244 

0.0440 

30.1139 

- 0.0246 

- 0.0695 

30.1987 

0.0030 

0.0244 

30.2836 

0.0082 

0. 0169 

30.3684 

0.0384 

0.0116 

30.4532 

0.0340 

- 0.0097 

30.5380 

0.0228 

0.0022 

30.6229 

0.0088 

0.0613 

30.7077 

0.0429 

- 0.0202 

30.7925 

0.0227 

0.0303 

30.8773 

0.0010 

0.0034 

30.9622 

- 0.0170 

0.0431 

3 i . 0470 

0.0416 

- 0.0445 

31.1318 

0.0289 

- 9.0104 

31.2167 

- 0.0319 

0.0054 

31.3015 

0.0149 

0.0131 

31.3863 

- 0.0026 

- 0.0054 

31.4711 

0.0081 

- 9.0199 

3 L . 5560 

- 0. 0069 

- 0.0321 

31.6408 

0.0073 

0.0244 

31.7256 

0.0493 

0.0110 

31.8105 

- 0. 0513 

- 0.0052 

31.8953 

0.0044 

- 0.0029 

31.9801 

0.0327 

0.0253 

32.0649 

- 0.0173 

0.0276 

32.1498 

0.0043 

- 0.0703 

32.2346 

- 0. 0053 

- 0 . OOCO 

32.3194 

0.0301 

0.0123 

32.4043 

0.0130 

- 0.0074 

32.4891 

0.0150 

- 0.0009 

32.5739 

0.0508 

0.0118 

32.6587 

0.0080 

0. 0333 

32.7436 

- 0.0221 

0.0249 

32.8284 

0 . C 592 

- 0.0363 

32.9132 

- 0.0183 

0.0234 

32.9980 

0.0148 

0.0148 

33.0829 

0.0219 

- 0.0376 

33.1677 

0.0549 

- 0.0027 

33.2525 

- 0.0344 

0.0020 

33.3374 

0.0115 

- 0.0002 

33.4222 

0.0392 

- 0.0086 

33.5079 

0.0156 

- 9.0226 

33.5918 

0.0173 

0.0266 

33.6767 

0 . C 009 

0.0295 

33.7615 

0.0395 

- 0.0034 

33.8463 

- 0 . C 196 

0.0522 

33.9312 

- 0.0073 

- 0.0178 
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g_ _ .079 S 3 + 11.0 S 2 + 60.0 S - 2.04 
S R D 2 


1 3_ = - 3.48 S 2 + 80.0 S + 33. I 
5 A d 2 

= ~ 0-168 S 3 - 10.37 S 2 - 54.64 S + 0.868 

H ' °2 

£_ = 1-35 S 2 - 12.8 S - 58.9 
S R 0 2 

= 204.0 S 2 + 237.0 S + 1883.0 
6 A d 2 


<£_ = 7.07 S 2 - 5.54 S 


where: 


D 2 = 0.998 S 4 + 6.70 S 3 + 15.8 S 2 + 53.9 S - 0.868 


= 0.998 (S + 5.6) (S - 0.0 16) (S + 0.564 + j 3.062HS + 0.564 - 


TABLE- III. Open Loop Lateral Transfer Functions 


i 3.062) 
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I .01 


U 0 = 198 ft/sec 
g = 32.2 ft/sec 2 

0 O = 0.0 
Y v = - 0. 167 

L p " - 5-54 

L p = I .35 
Lg = - 7.06 
N p = - 0. 156 


Ng = 9.27 

Y* = 0.0 

°A 

Y* = 0.079 

°R 

Is = 204.0 
°A 

l 6 = L95 

°R 

Hs = - 2.99 
°A 

N. = -10.7 
6 R 


TABLE IV. Lateral Equations of Motion Parameters 


U Q = 198 ft/sec 
g = 32.2 ft/sec 2 

X u = - 0.0466 

X w = 

X w = 0.0793 

Xq = 0-0 

Z u = - 0.325 
Z; = - 0.0154 
Z w = - 3.045 

Z q “ - 8 - 12 

M u = 0.00079 

M; = - 0.0106 


M w = ~ 0.175 
M q = ~ 5.56 

X$ s = 0.0 

X 5f _ = - 8.84 

X a = 15.7 

. y 

Z 6s = 89.3 

Z($ F = -283.1 

2 n = - 3.05 
9 

\ = 61.2 
% = -17.8 
M„ = 0.156 

a g 


TABLE V, Longitudinal Equations of Motion Parameters 


a 

- 1.433 S 3 

- 25.2 S 2 - 3.814 S - 1 .082 

<*F = 


D l 

a _ 

0.452 S 3 + 

61.29 S 2 + 9.49 S + 3.354 

6 S 


D l 

a 

- 3.045 S 3 

- 50.62 S 2 - 7.876 S - 2. 131 

a g = 


D l 

e _ 

- 15.03 S 2 

- 6.851 S - 1.805 

% 


D l 

e_ _ 

61.22 S 2 + 

180.3 S + 28. i 

6 s 


D l 

Q_ = 

- 28.78 S 2 

- 4.364 S 

a 

g 

D | 

1 


where: 

D, = 1.015 S 4 + 10.85 S 3 + 51.80 S 2 + 8.075 S + 2.131 

= 1.015 (S + 0.076 + j 0 . 1916) (S + 0.076 - j 0.1916) 

x (S + 5.268 + j 4.65) (S + 5.268 - j 4.65) 


TABLE VI , Open Loop Longitudinal Transfer Functions 
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FIGURE I. TIME HISTORY OF VERTICAL GUST USED TO EXCITE 
AIRCRAFT 







FIGURE 2. TIME HISTORY OF LATERAL GUST USED TO 
EXCITE AIRCRAFT 






2 


3 


FIGURE 


FREQUENCY - LOG a) 

APPROXIMATE ONE -DIMENSIONAL POWER 
SPECTRAL DENSITY OF LATERAL GUST 




OPEN LOOP 



RESPONSE 
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FIGURE 5 A. ANGLE OF ATTACK RESPONSE TO 
VERTICAL GUST INPUT 


TIME (SECONDS) 


CLOSED LOOP 



TIME (SECONDS) 




OPEN LOOP 



FIGURE 5C. PITCH ANGLE RESPONSE TO 
36 VERTICAL GUST INPUT 


TIME (SECONDS) 


CLOSED LOOP 



TIME (SECONDS) 








CLOSED LOOP 



FIGURE 6B. AERODYNAMIC SIDESLIP ANGLE 

RESPONSE TO LATERAL GUST INPUT 
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TIME (SECONDS) 







FIGURE 6D. ROLL ANGLE RESPONSE TO 
LATERAL GUST INPUT 


TIME (SECONDS) 
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FIGURE 7 UNITY FEEDBACK REPRESENTATION 
OF LATERAL CONTROL SYSTEM 





FIGURE 8 MODIFIED AIRCRAFT 
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FIGURE 9. UNITY FEEDBACK REPRESENTATION OF 
LONGITUDINAL CONTROL SYSTEM 
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